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Abstract: 

Fast  responding  FP-07  thermistors  have  been  incorporated  on  profiling  EM-APEX  floats  to 
measure  microscale  ocean  temperature  fluctuations  produced  by  turbulence.  In  this 
implementation,  the  FP-07  thennistor  generates  an  electrical  signal  corresponding  to  ocean 
temperature  fluctuations,  which  is  conditioned  by  an  analog  circuit  board,  and  digitized  and 
recorded  to  a  custom  data  acquisition  and  storage  board.  The  raw  and  processed  temperature 
observations  are  stored  on  a  microSD  card.  Results  from  eight  microstructure  EM-APEX  floats 
deployed  in  the  Sargasso  Sea  are  presented  here.  The  slow  profiling  speed  of  EM-APEX  floats 
enables  them  to  capture  the  higher  wavenumber  regime  of  microscale  temperature  variation.  The 
quality  of  temperature  variance  dissipation  rates  %  estimated  from  microstructure  EM-APEX 
floats  is  verified  by  their  agreement  with  the  Batchelor  spectrum  and  by  the  close  inter-float 
agreement  of  temporal  and  vertical  variations  measured  by  multiple  floats.  Estimates  of  %  from 
the  profiling  floats  exhibit  a  lognormal  distribution  as  expected  for  statistically  homogeneous 
isotropic  turbulence.  Turbulence  measurements  derived  from  FP-07  sensors  on  autonomous 
profiling  floats  are  of  comparable  quality  to  those  on  conventional  free-fall  microstructure 
profilers. 
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1.  Introduction 

Quantifying  oceanic  turbulence,  diffusion,  and  dissipation  gives  insight  to  important  oceanic 
processes  and  improves  parameterization  schemes  in  numerical  models.  The  intermittency  of 
turbulence  and  the  difficulty  obtaining  adequate  realizations  of  turbulence  observations  limit  the 
reliability  of  turbulence  estimations.  One  sensor,  the  FP-07  thennistor,  is  capable  of  rapid 
temperature  sampling  at  microscales,  and  has  been  used  extensively  on  free-fall  vertical  profilers 
and  moorings  ( Gregg  and  Meagher,  1980;  Mourn  and  Nash,  2009).  The  temperature  variance 
dissipation  rate  x  can  be  computed  from  FP-07  microscale  temperature  measurements.  The 
development  of  autonomous  vertical  profilers  and  the  ability  to  deploy  several  in  an  array  allows 
collection  of  many  profiles  at  multiple  sites  simultaneously.  Because  of  the  dynamical  link 
between  turbulence  and  small-scale  vertical  shear  and  strain,  autonomous  profiling  floats  with 
FP-07  thennistors  as  well  as  other  sensors  to  observe  conventional  water  properties  are  ideal  to 
better  measure  and  understand  small-scale  and  microscale  processes  and  improve  physics-based 
turbulence  parameterization  schemes.  Here,  we  describe  our  approach  to  implement  turbulence 
sensing  on  the  EM-APEX  float  and  present  results  from  deployments  of  multiple  floats. 


2.  Material  and  methods 

2.1  EM-APEX  floats 

The  EM-APEX  float  combines  the  standard  Teledyne  Webb  Research  Corp.  APEX  profiling 
float  with  an  APL-UW  subsystem  that  measures  the  motion-induced  electric  fields  generated  by 
the  ocean  currents  moving  through  the  vertical  component  of  the  earth’s  magnetic  field  (Fig.  1) 
{Sanford  et  al.,  2005).  The  APEX  float  changes  its  buoyancy  to  enable  it  to  profile  the  ocean 
with  a  maximum  depth  range  of  2000  m.  Temperature  and  salinity  measurements  are  obtained 
from  a  Sea-Bird  Electronics  SBE41CP  CTD,  and  taken  every  10-25  s.  When  on  the  sea  surface, 
the  float’s  position  is  determined  by  GPS,  and  the  position  and  accumulated  profile  data  (except 
for  the  raw  FP-07  data)  are  transmitted  over  the  Iridium  global  satellite  phone  system. 

The  velocity  sensor  operates  on  the  same  principles  of  motional  induction  applied  on  the 
Absolute  Velocity  Profiler  {Sanford  et  al.,  1985)  and  the  Expendable  Current  Profiler  {Sanford  et 
al.,  1982).  The  electric  field  sensing  electrodes  are  located  on  the  top  end  of  the  floats.  Other 
necessary  measurements  are  magnetic  compass  heading  and  instrument  tilt.  The  float  descends 
and  ascends  at  ~0. 1 5  m  s~ 1  and  rotates  at  a  period  of  ~12  s.  The  motional-induced  electric  field  is 
determined  by  a  sinusoidal  fit  to  the  measured  voltages  using  the  basis  functions  from  the 
horizontal  components  of  the  magnetometer.  The  fit  is  made  over  a  50-s-long  segment  of  data, 
and  the  averaging  window  is  moved  25  s  between  successive  fits.  These  sinusoidal  fits  and  the 
nns  (root  mean  square)  residuals  are  transferred  to  the  APEX  float  controller  for  storage  and 
later  transmission  over  Iridium.  The  fits  provide  an  estimate  of  the  horizontal  current,  and  the 
rms  residuals  provide  an  estimate  of  the  velocity  noise  level.  In  the  upper  ocean,  the  rms 
residuals  can  also  be  used  to  estimate  properties  of  surface  gravity  waves,  e.g.,  amplitude,  period, 
and  wavelength.  The  horizontal  current  uncertainty  is  about  0.015  ms1. 


2.2  Turbulence  sensor  implementation 

A  fast-response  temperature  sensor,  the  FP-07  thermistor,  and  electronic  signal  conditioning 
board  are  provided  by  Rockland  Scientific,  Inc.  (RSI),  and  combined  with  a  custom  digital  data 
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acquisition  board.  The  complete  unit  is  installed  in  EM -APEX  floats  (Fig.  1).  In  addition  to  the 
thermistors  installed  in  metal  shafts  or  stings,  the  implementation  requires  special  end  caps  for 
the  EM- APEX  float. 

The  arrangement  of  the  major  components  within  the  EM- APEX  float  is  outlined  in  a  block 
diagram  (Fig.  2).  Temperature  fluctuations  cause  the  FP-07  thermistor  to  change  its  resistance. 
Electric  current  passing  through  the  FP-07  thennistor  produces  voltage  fluctuations  caused  by 
these  resistance  fluctuations.  The  voltages  are  amplified  and  conditioned  on  the  RSI  pASTP-LP 
board  and  passed  to  the  APL-UW  board.  The  APL-UW  board  (MLADC)  digitizes  the  amplified 
voltage  signals  to  24-bit  resolution.  The  board  also  contains  a  low-power  microprocessor  that 
buffers  and  stores  the  digital  data  onto  a  microSD  card.  In  the  latest  version  of  the  microstructure 
package,  temperature  spectra  are  computed  and  band-averaged  spectra  are  transmitted  over 
Iridium.  The  turbulence  circuitry  is  powered  and  controlled  by  the  APEX  controller  APF-9 
board. 


2.3  Estimates  of  %,  turbulent  temperature  diffusivity  Kj  and  turbulent  kinetic  energy  dissipation 
rate  s 


The  high-frequency  FP-07  sensor  provides  measurements  of  turbulent  temperature 
fluctuations  that  can  be  used  to  compute  the  temperature  variance  dissipation  rate  %,  turbulent 
temperature  diffusivity  Kj,  and  turbulent  kinetic  energy  dissipation  rate  s.  The  temperature 
variance  dissipation  rate  for  isotropic  turbulence  is  defined  as 


*  =  ^  <©2  +  (|)2  +  ©2>  =  6  Dr  <©2>  =  6  Dt  /”  *,<*,)**, 


where  DT  is  the  molecular  temperature  diffusivity  (-1.4  x  10  7  nr  s  1 ),  <pgzT(.kz)  the  vertical 

wavenumber  kz  spectrum  of  the  vertical  gradient  of  temperature  and  k0  the  lowest 

wavenumber  of  significant  turbulence.  The  FP-07  measures  vertical  gradients  of  temperature 
well  into  the  turbulent  dissipation  wavenumber  regime. 


Three  transfer  functions  are  applied  to  FP-07  sensor  measurements  to  derive  the  temperature 
gradient  spectrum;  they  are  the  analog-to-digital  conversion  Hadc,  amplification  of  the  high- 
frequency  component  HPreamp,  and  the  response  of  the  sensor  to  temperature  fluctuations  HFP07, 
where  H  represents  the  transfer  function.  The  analytical  fonns  of  the  transfer  functions  Hadc  and 
HPreamp  are  provided  by  manufacturers  Cirrus  Logic  and  Rockland  Scientific,  respectively. 

The  transfer  function  HFP07  is  not  known  exactly  and  different  findings  have  been  reported 
(Nash  et  al.  1999;  Sommer  et  al.  2013).  We  examine  the  effects  of  three  double-pole  transfer 
functions  proposed  by  Vachon  and  Lueck  ( 1 984),  Gregg  and  Meagher  ( 1 980),  and  Sommer  et  al. 
(2013)  on  our  estimates  of  %.  The  squared  double-pole  transfer  function  (H|P07)  is  expressed  as 


H F P 07  (.fY 


1 

[1  +  (27t/t)2]2 


where/is  the  frequency  in  Hz  and  x  is  a  constant  characterizing  the  time  in  seconds  of  FP-07 
sensor  response  to  ambient  temperature  fluctuation.  The  time  constant  is  suggested  to  vary  with 
the  flow  speed  relative  to  the  sensor  expressed  as, 
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where  r0  is  a  reference  time  constant,  and  U  the  flow  speed  relative  to  the  sensor  in  m  s  1  and  Uo 
=  1  m  s.  The  reference  time  r0  is  0.0064  s,  0.005  s,  and  0.010  s,  and  the  speed-dependence 
power  law  (s)  is  -0.5,  -0.32,  and  0  for  the  transfer  functions  suggested  by  Vachon  and Lueck 
(1984),  Gregg  and  Meagher  (1980),  and  Sommer  et  al.  (2013),  respectively.  At  the  typical  EM- 
APEX  float  vertical  profiling  speed  of  0.15  m  s”1,  the  time  constant  x  is  0.0155  s,  0.0092  s,  and 
0.010  s  for  the  Vachon  and  Lueck  (1984),  Gregg  and  Meagher  (1980),  and  Sommer  et  al.  (2013) 
transfer  functions,  respectively.  The  latter  two  transfer  functions  are  in  very  close  agreement. 

FP-07  sensors  are  mounted  on  the  top  end  of  EM- APEX  floats,  so  measurements  are  taken 
only  during  ascending  profiles  to  avoid  potential  wake  effects.  The  FP-07  samples  temperature  at 
120  Hz.  Estimates  of  x  were  computed  from  4096  temperature  data  points  over  ~34-s  duration, 
corresponding  to  ~5-m  vertical  bin  for  a  typical  EM- APEX  float  ascending  speed  of  0.15  m  s”1. 

Computing  x  from  one  timeseries  of  34-s  duration  begins  with  the  voltage  measured  by  the 
FP-07  converted  to  temperature,  then  calibrated  with  the  standard  SBE41CP  CTD  sensor  on  the 
float  (Fig.  3).  The  frequency  spectrum  of  temperature  is  computed  and  converted  to  the  vertical 
wavenumber  spectrum  assuming  a  frozen  hypothesis  using  the  average  ascending  float  speed 
within  the  34-s  time  interval  (Fig.  3b).  The  vertical  temperature  gradient  spectrum  is  computed 
as  vertical  wavenumber  squared  times  the  vertical  wavenumber  temperature  spectrum  (Fig.  3c). 
The  vertical  temperature  gradient  spectrum  is  then  corrected  by  the  above-mentioned  three 
transfer  functions.  The  Hadc  squared  is  white  at  low  wavenumbers  and  rolls  off  beyond  ~5  x  10“ 
m”1.  The  Hpreamp  squared  is  white  at  the  lowest  wavenumbers,  has  a  +2  spectral  slope  beyond  ~  10 
in”1,  and  rolls  off  at  wavenumbers  beyond  that  shown  in  Fig.  3e.  Transfer  functions  proposed  by 
Sommer  et  al  ( 1 984)  and  Gregg  and  Meagher  ( 1 980)  are  in  close  agreement,  but  that  proposed 
by  Vachon  and  Lueck  (1 9 84)  rolls  off  at  a  frequency  about  30%  less  (Fig.  3d). 

The  vertical  temperature  gradient  spectra  corrected  by  the  three  HFP07  transfer  functions  are  in 
close  agreement  (Fig.  3e-f).  Because  of  the  low  profiling  speed,  the  roll-off  spectral  shape  of  the 
temperature  gradient  spectrum  is  fully  resolved,  implying  that  most  turbulent  temperature 
gradient  variance  is  captured.  The  noise  spectrum  shows  a  spectral  shape  of  kz2  at  wavenumber 
greater  than  10  m  . 

Estimates  of  x  are  computed  by  integrating  the  temperature  gradient  spectrum  from  kz  =  12 
m”1  to  an  upper  bound  kz  beyond  which  the  observed  temperature  gradient  spectrum  is  smaller 
than  the  noise  spectrum.  The  lower  bound  wavenumber  12  m”1  is  detennined  after  examining 
observed  spectra.  Estimated  x  in  this  example  is  4.9  x  10  K  s  using  the  FP-07  transfer 
function  of  Vachon  and  Lueck  (1984)  and  3.6  x  10  9  K2  s”1  using  those  of  Gregg  and  Meagher 
(1980)  and  Sommer  et  al.  (2013).  Estimates  of  x  derived  from  different  transfer  functions  never 
differ  by  more  than  a  factor  of  two.  Note  that  because  the  full  turbulent  temperature  gradient 
spectrum  is  resolved,  estimates  of  x  are  computed  by  integrating  the  observed  temperature 
gradient  spectrum  over  the  turbulence  range  without  the  assumption  of  a  model  temperature 
gradient  spectrum,  as  proposed  by  Batchelor  (1959)  or  Kraichnan  (1967).  Observed  temperature 
gradient  spectra  are  in  good  agreement  with  the  Batchelor  (1959)  temperature  spectrum. 

The  turbulence  diffusivity  of  temperature  KT  can  be  computed  following  Osborn  and  Cox 
(1972) as 
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where  (— )  is  the  background  temperature  gradient.  Assuming  that  the  turbulence  diffusivity  of 

temperature  KT  equals  the  turbulent  diffusivity  of  buoyancy  (e.g.,  Mourn  and  Nash,  2009;  Alford 
and  Pinkel,  2000),  defined  as  Kp  =  Te/N2,  where  £  is  the  turbulent  kinetic  energy  dissipation 
rate,  N  the  buoyancy  frequency,  and  T  the  mixing  efficiency  assumed  as  the  nominal  value  of  0.2, 

N2y 

we  can  estimate  £  as  £  =  — -tf—  ■ 


3.  Results 

3.1  Experimental  observations 

Twenty  EM-APEX  floats  were  deployed  in  the  Sargasso  Sea  southeast  of  Cape  Hatteras  in 
2011  as  part  of  the  ONR- funded  Scalable  Lateral  Mixing  and  Coherent  Turbulence  Departmental 
Research  Initiative  (Fig.  4)  ( Shcherbina  et  al.,  2015).  All  20  EM-APEX  floats  measured 
temperature,  salinity,  pressure,  and  horizontal  velocity  components  between  the  surface  and 
typically  150-m  depth,  with  some  deeper  profiles  to  ~250  m.  Although  the  floats  were  at 
different  locations,  with  separations  up  to  10  km,  the  units  were  programmed  to  rise  nearly  in 
synchronism.  Extensive  effort  was  made  to  synchronize  the  profiling  floats.  Simultaneous 
profiles  eliminate  differences  caused  by  temporal  versus  spatial  differences.  Ten  of  the  EM- 
APEX  floats  were  equipped  with  dual  FP-07  sensors.  One  of  these  % -EM-APEX  floats  failed 
during  the  first  deployment  and  was  not  used  again  in  the  experiment,  and  another  was  lost  in  the 
second  deployment.  Therefore,  eight  ^-EM-APEX  floats  provided  turbulence  measurements. 

Floats  were  deployed  in  three  settings  (Fig.  4):  Deployment  (Dl)  was  in  a  “weak  energy” 
region  during  3-10  June  2011  (4988  profiles),  Deployment  2  (D2)  was  in  a  “strong  front” 
during  13  -  16  June  2011  (2546  profiles),  and  Deployment  3  (D3)  was  conducted  in  a  “strong 
current”  during  17-20  June  2011  (2032  profiles)  (Fig.  4).  Floats  were  intended  to  be  deployed 
in  three  concentric  circles  of  0.5  m,  1  km,  and  2  km  radii  in  the  first  two  deployments,  and  of  1,  2, 
and  4  km  radii  in  the  third  deployment,  with  roughly  six  floats  on  each  circle.  However,  because 
of  the  oceanic  horizontal  shear,  float  arrays  were  not  deployed  in  a  perfectly  concentric 
configuration. 

During  the  7-d  mission  of  the  first  deployment,  20  EM-APEX  floats  were  advected  southeast 
about  20  km  (Fig.  4).  There  was  little  dispersion  of  the  float  array.  During  the  3-d  mission  of  the 
second  deployment,  the  float  array  was  strongly  sheared  in  a  front,  while  moved  northeast  for 
-30  km.  During  the  3-d  mission  of  the  third  deployment,  the  float  array  was  advected  northeast 
for  100  km. 


3.2  Turbulence  properties 

The  methods  for  computing  y,  KT,  and  s  are  discussed  in  section  2.3.  Because  KT  and  s  are 
derived  quantities  from  measurements  of  y  and  our  purpose  is  to  report  the  quality  of  turbulence 
measurements  derived  from  FP-07  sensors  on  the  EM-APEX  floats,  we  will  discuss  only  the 
quality  of  %  in  the  following  analysis.  The  quality  of  %  will  be  assessed  by  three  independent 
methods:  (1)  comparing  the  observed  temperature  gradient  spectrum  with  the  canonical  fonn  of 
the  turbulent  temperature  gradient  spectrum  ( Batchelor ,  1959),  (2)  comparing  estimates  of  % 
obtained  from  different  floats,  and  (3)  comparing  the  statistics  of  %  estimates  with  the 
expectation  for  equilibrium  turbulence. 
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In  the  following  analysis,  estimates  of  %  computed  using  the  transfer  function  of  Gregg  and 
Meagher  (1984)  will  be  discussed.  The  ratio  of  estimates  of  %  computed  from  two  FP-07  sensors 
has  a  mean  ratio  of  1.06,  and  the  95%  confidence  interval  between  0.50  and  1.88.  Estimates 
X  with  the  ratios  between  those  obtained  from  two  FP-07  sensors  less  a  factor  of  two  are 
averaged  and  discussed  in  the  following  analysis. 

Observed  temperature  gradient  spectra  and  their  variance  preserving  spectra  are  compared 
with  those  of  Batchelor  ( 1959)  at  %  of  10  10,  10"9,  10  8,  10  7,  and  10  6  K2  s_1  (Fig.  5).  Note  that 
each  observed  temperature  gradient  spectrum  is  associated  with  distinct  values  of  %  and  s  .  The 
Batchelor  spectrum  is  computed  for  each  corresponding  observed  spectrum  with  its  associated 
X  and  s.  These  Batchelor  spectra  are  averaged  at  different  bins  of  %  and  compared  with  the 
averaged  observed  spectra.  Observed  spectra  agree  well  with  Batchelor  spectra  near  the  roll-off 
wavenumbers,  where  most  temperature  gradient  variance  resides.  They  differ  at  low 
wavenumbers,  where  the  temperature  gradient  variance  is  small  and  therefore  does  not  contribute 
significantly  to  the  estimate  of  %.  The  noise  level  of  our  %  estimates  is  about  10  K“  s  ,  where 
the  peak  of  the  variance  is  closer  to  the  noise.  The  noise  is  likely  due  to  the  amplified  electronic 
noise  in  the  float. 

Vertical  profiles  of  %  estimated  from  eight  floats,  averaged  over  each  of  three  deployments, 
show  similar  vertical  structures  with  their  magnitudes  typically  within  a  factor  of  three  from  the 
means  from  all  floats  (Fig.  6).  Estimates  of  %  obtained  from  different  floats  at  different  depths 
are  strongly  correlated  with  x  averaged  over  all  floats  with  a  correlation  coefficient  R  =  0.94, 
much  higher  than  the  95%  significance  level  of  0.07  (Fig.  7).  The  variation  is  less  than  a  factor 
of  two  for  x  greater  than  3x10  K  s  and  a  factor  of  six  for  smaller  x,  where  the  noise 
becomes  significant  (Fig.  5).  The  temperature  signal  is  contaminated  by  electronic  noise  in  the 
float.  In  a  later  version  of  the  x-EM-APEX  system,  the  preamps  are  relocated  outside  the  endcap 
to  minimize  the  contamination  by  electronic  noise  from  within  the  float. 

The  temporal  variation  of  depth-averaged  x  observed  by  different  floats  is  similar  most  of  time 
(Fig.  8).  The  greatest  variability  of  depth-averaged  x  ,  observed  during  15-16  June,  is 
presumably  associated  with  a  real  oceanic  signal  of  turbulence  variability  when  the  float  array 
encountered  a  front  ( Shcherbina  et  al.,  2015).  At  other  times,  the  variability  of  the  depth 
averaged  %  of  different  floats  is  about  a  factor  of  two. 

For  statistically  homogeneous  isotropic  turbulence,  we  expect  s  and  x  to  have  a  lognonnal 
distribution  (Yaglom,  1960;  Gregg  et  al.,  1993;  Yamazaki,  1990).  The  probability  density 
function  of  x  estimates  computed  from  all  floats  shows  a  nearly  lognonnal  distribution,  with  a 
mean  logio(x)  of -8.4  and  a  standard  deviation  of -0.9  (Fig.  9).  The  Quantile-Quantile  (Q-Q) 
plot  shows  the  quantile  of  observed  logio(x)  vs.  the  standardized  nonnal  distribution  ( Chambers 
et  al.  1983).  The  agreement  is  good  for  logio(x)  between  -9.5  and  -6.8,  within  ~1.2  standard 
deviation  from  the  mean.  Yamazaki  and  Lueck  (1990)  report  that  when  turbulence  is  not 
statistically  homogeneous,  the  observed  turbulent  kinetic  energy  dissipation  rate  £  might  not 
exhibit  a  lognormal  distribution.  We  suspect  that  the  deviation  of  lognonnal  distribution  at  low 
X  is  due  to  the  measurement  noise,  and  the  deviation  at  high  %  may  be  due  to  a  violation  of  the 
assumptions  of  statistical  homogeneity. 

4.  Discussion  and  conclusions 
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Our  assessment  of  the  quality  of  %  estimates  is  based  on  three  independent  methods. 
Estimates  are  reliable  at  values  greater  than  order  of  10  K~s  .  Performance  is  at  least 
comparable  to  that  achieved  in  conventional  microstructure  free-fall  profilers,  such  as 
X  estimates  from  Advanced  Microstructure  Profiler  observations  ( Lien  et  al.,  1995).  Most 
importantly,  because  of  the  float’s  slow  profiling  speed,  the  roll-off  regime  of  the  turbulence 
gradient  spectrum  is  captured.  In  comparison,  free-fall  microstructure  profilers  typically  fall  at  a 
speed  of  0.5-1  ms"1,  3-7  times  faster  than  the  float’s  profiling  speed,  in  order  to  capture 
turbulence  shear  variance  measured  from  air-foil  shear  probes.  Unfortunately,  at  this  high  speed, 
the  roll-off  wavenumber  of  the  temperature  gradient  spectrum  is  often  not  captured  by  the  FP-07 
sensor. 

The  installation  of  FP-07  sensors  on  the  eight  EM -APEX  floats  used  in  the  LatMix  swann 
was  able  to  observe  the  small-scale  temperature  fluctuations  associated  with  ocean  turbulence. 
The  observations  conform  to  the  expectations  of  the  lognonnal  distribution  of  homogeneous 
turbulence  at  the  intermediate  levels  of  turbulence.  Significant  differences  in  turbulence 
temperature  dissipate  rates  were  calculated  from  floats  across  a  front.  This  is  expected  in  a  region 
of  spatial  gradients  in  oceanic  properties. 

The  performance  of  the  %  system  was  adequate  for  this  first  application  on  the  EM-APEX 
floats.  There  were  some  perfonnance  limitations.  A  concurrent  installation  of  the  RSI  FP-07  and 
preamp/conditioning  board  on  Seagliders  by  our  colleagues  Luc  Rainville  and  Jason  Gobat  has  a 
much  lower  noise  level.  This  suggests  that  our  design  is  contaminated  by  the  amplified  EM  noise 
inside  the  float.  The  stainless  steel  metal  stings  also  contributed  corrosion  noise  that  was 
observed  on  the  Ag-AgCl  electrodes  near  the  end  cap.  These  led  to  noise  in  the  float’s  velocity 
measurements.  Both  of  these  effects  have  been  reduced  in  the  latest  design.  The  thennistor 
preamps  are  now  within  the  end  cap,  so  that  electronic  noise  within  the  float  is  not  amplified  and 
its  contamination  is  minimized.  Also,  the  stainless  steel  stings  have  been  replaced  with  an  all 
plastic  holder  for  thennistor  probes  and  preamps. 

The  need  to  recover  the  floats  to  extract  data  on  the  microSB  card  prevented  the  use  of  this 
design  in  other  experiments.  The  latest  firmware  has  been  changed  to  provide  log  averaged 
spectral  values.  The  amount  of  data  to  be  transmitted  over  Iridium  is  expected  to  be  about  the 
same  order  as  that  for  the  CTD  and  velocity  data. 
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Figure  1 :  Photos  of  EM-APEX  float  with  electrode  collar  and  rotation  vanes  (left)  with  dual  fast- 
response  temperature  sensors  FP-07  inside  a  guard  (light  green)  and  (right)  a  closer  view  of  the 
FP-07  sensors,  visible  on  black  tips  of  metal  shafts.  Other  components  on  the  top  of  the  float  are 
the  SBE41CP  CTD  in  the  center  and  the  GPS/Iridium  antenna. 
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Figure  2:  Block  diagram  of  the  /-EM-APEX  float.  The  left  side  of  the  dashed  line  shows  the 
standard  configuration  for  an  EM-APEX  float  (grey).  The  APF9i  APEX  float  controller  runs  the 
Sea-Bird  CTD  while  using  the  air  pump,  air  valve,  and  piston  motor  to  control  buoyancy.  At  the 
surface  the  APF9i  uses  a  single  antenna  with  an  RF  switch  for  GPS  position  information  and 
Iridium  communications.  The  right  side  of  the  dashed  line  shows  the  components  added  to  make 
an  '/-EM-APEX  float  (light  green),  including  dual  RSI  FP-07  fast  temperature  probes  (dark 
green),  differential  RSI  pre-amps  with  a  DC  gain=x6  feeding  into  a  24-bit  sigma-delta  analog-to- 
digital  converter  and  a  dedicated  /-EM-APEX  controller.  The  APF9i  controls  power  to  the  / 
system  using  the  persistent  power  interface  board.  The  /  controller  on  the  MLADC  board  uses 
the  CTD  sniffer  circuit  board  to  intercept  the  data  stream  from  the  CTD  to  the  APF-9i  to  gather 
rise  and  fall  rate  data.  The  resulting  /-EM-APEX  data  are  stored  on  a  MicroSD  card.  Reduced 
data  are  transmitted  over  the  satellite  link  when  on  the  surface,  but  the  raw  data  are  available 
only  by  removal  of  the  storage  card. 
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Figure  3:  Examples  of  computing  temperature  variance  dissipation  rate  %  from  FP-07 
measurements:  (a)  the  time  series  of  35  s  of  temperature  from  one  FP-07  sensor,  (b)  the  observed 
raw  temperature  spectrum,  (c)  the  observed  raw  spectrum  of  vertical  gradient  of  temperature,  (d) 
squared  transfer  functions  of  pre-amplifier  (HPreamp,  green  curve  labeled  as  Preamp),  analog-to- 
digital  converter  (Hadc,  ADC)  (black  curve  labeled  as  ADC),  FP-07  sensor  (FIFP07,  thick  blue 
curve  labeled  as  FP07vl  and  thin  blue  curve  labeled  as  FP07gm),  and  the  product  of  the  above 
three  transfer  functions  (thick  red  curve  labeled  as  Total vl  and  thin  red  curve  labeled  as  TotalGM), 
(e)  observed  temperature  gradient  spectrum  corrected  for  transfer  functions  (thick  and  thin  black 
curves  labeled  as  Obs)  and  Batchelor  temperature  gradient  spectrum  (thick  and  thin  red  curve 
labeled  as  Bat),  and  (f)  the  variance  preserving  form  of  spectra  shown  in  panel  (e).  The  subscript 
‘VL’  represents  the  use  of  FP-07  transfer  functions  proposed  by  Vachon  and  Lueck  (1984).  The 
subscript  ‘GM’  represents  the  use  of  the  transfer  function  suggested  by  Gregg  and  Meagher 
(1980),  which  is  nearly  identical  to  that  of  Sommer  et  al.  (2013).  Note  that  Batchelor  spectra 
shown  in  panels  (e)  and  (f)  are  computed  using  estimates  of  x  and  £  obtained  from  observed 
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spectra.  Thick  and  thin  curves  in  panels  (e)  and  (f)  are  computed  using  FP-07  transfer  functions 
of  Vachon  and  Lueck  ( 1 984)  and  Gregg  and  Meagher  (1980),  respectively. 
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Figure  4:  Float  trajectories  of  three  deployments  during  the  LatMix  experiment.  Blue  curves 
represent  float  trajectories  on  6-10  June  (Dl),  black  curves  on  13-16  June  (D2),  and  pink  curves 
on  17-20  June  (D3).  Background  color  shows  the  AVFIRR  image  of  SST  on  4  June. 
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Fig.  5:  (a)  Comparison  among  observed  temperature  gradient  spectra  (dashed)  with  Batchelor 
spectra  (solid)  for  x  of  10  10  -  10  6  K2  s_1.  (b)  Similar  to  (a)  but  in  variance  preserving  format,  (b) 
Gray  shadings  represent  95%  confidence  interval  of  observed  spectra,  and  thin  curves  represent 
95%  confidence  interval  of  Batchelor  spectra. 
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Vertical  Profiles  of  Temporally  Averaged  x 


Figure  6.  Vertical  profiles  of  x  temporally  averaged  over  three  deployments  of  different  floats 
(a-c)  and  averaged  over  all  floats  (d).  In  panels  (a)-(c)  different  thin  colored  curves  represent 
results  from  different  floats,  and  the  thick  curves  represent  averages  over  all  floats. 
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Fig.  7:  Scatter  plot  of  %  temporally  averaged  over  each  deployment  period  for  different  floats  vs. 
X  temporally  averaged  over  each  deployment  period  and  over  all  floats.  Colors  represent 
different  floats. 
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Time  Series  of  Depth  Averaged  % 
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Fig.  8:  (a)  Time  series  of  depth  averaged  %  from  all  floats,  (b)  six -hourly  averaged  values. 
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Fig.  9:  (a)  Probability  density  function  of  the  ensemble  of  logio(x)  observations  from  all  floats 
(black  curve)  compared  with  that  expected  from  a  theoretical  normal  distribution  of  the  observed 
mean  (-8.3)  and  standard  deviation  (0.9)  (red  curve),  (b)  Q-Q  plot  of  logio(x)  vs.  that  of 
standardized  nonnal  distribution.  The  blue  curve  shows  observed  Q-Q  structure.  Red  curves 
represent  the  95%  confidence  interval  assuming  a  normal  distribution.  The  black  dashed  curve 
represents  the  theoretical  normal  distribution. 
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